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Abstract 

The product of the DKC1 gene, dyskerin, is required for both ribosome biogenesis and telomerase complex stabilization. 
Targeting these cellular processes has been explored for the development of drugs to selectively or preferentially kill cancer 
cells. Presently, intense research is conducted involving the identification of new biological targets whose modulation may 
simultaneously interfere with multiple cellular functions that are known to be hyper-activated by neoplastic 
transformations. Here, we report, for the first time, the computational identification of small molecules able to inhibit 
dyskerin catalytic activity. Different in silico techniques were applied to select compounds and analyze the binding modes 
and the interaction patterns of ligands in the human dyskerin catalytic site. We also describe a newly developed and 
optimized fast real-time PCR assay that was used to detect dyskerin pseudouridylation activity in vitro. The identification of 
new dyskerin inhibitors constitutes the first proof of principle that the pseudouridylation activity can be modulated by 
means of small molecule agents. Therefore, the presented results, obtained through the usage of computational tools and 
experimental validation, indicate an alternative therapeutic strategy to target ribosome biogenesis pathway. 
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Introduction 

One basic prerequisite for the development of antineoplastic 
therapeutics is represented by the identification of cellular 
processes that are selectively altered in cancer cells and that could 
be modulated by pharmacological actions on specific biological 
targets. 

Among a series of cellular processes, both ribosome production 
and telomerase functions are known to be hyper-activated by 
neoplastic transformation. On one side, the rate of ribosome 
biogenesis regulates cellular growth and proliferation, and cancer 
cells carry over an increased production of ribosomes to sustain the 
protein synthesis necessary for unbridled cell growth [1,2]. On the 
other hand, the reactivation of telomerase, allowing the mainte- 
nance of chromosome ends during cell proliferation, is a 
characteristic of about 85-90% of primary tumors. Although it 
is not detectable in most somatic cells, with the exception of some 
adult pluripotent stem cells, proliferative cells of renewal tissues, 
and male germline cells [3,4]. Therefore the specific targeting of 
each of these two cellular processes has been explored for the 
development of drugs in order to selectively or preferentially kill 
cancer cells [5-8]. 

The product of DKC1 gene, dyskerin, is necessary for both 
processes of ribosome biogenesis and telomerase complex stabi- 



lization [9]. Indeed, dyskerin mediates the site specific uridine 
conversion to pseudouridine in rRNA and snRNA. Uridine 
modification in rRNA represents an early and crucial step of 
rRNA processing affecting the rate and the efficiency of ribosome 
production [10,1 1]. In addition, dyskerin also binds the telomerase 
RNA component (TERC), stabilizing the telomerase enzymatic 
complex and the mutations of the DKC 1 gene at specific sites or 
the reduction of its expression strongly reduces the levels of TERC 
and the activity of telomerase [9]. In the rare multisystemic 
syndrome X-linked dyskeratosis congenital and in a subset of 
human tumors arising in the general population, dyskerin has been 
proposed to act as a tumor suppressor [9,10,12]. In contrast, 
dyskerin is overexpressed in a number of human cancer types and 
high levels of dyskerin expression in tumors are associated with an 
aggressive clinical behavior in various tumor types including breast 
[10], prostate [11], head and neck [13], colon [14], and 
hepatocellular carcinomas [15]. These contrasting observations 
may be explained considering that in some cases the partial lack of 
dyskerin function could in the long term promote peculiar 
neoplastic features, while in a distinct subset of aggressive tumors 
the need to support the increased production of ribosomes and the 
increased demand for telomerase function characterizing actively 
growing tumor cells requires dyskerin overexpression. 
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This is consistent with the need to support the increased 
production of ribosomes and the increased demand for telomerase 
function characterizing actively growing tumor cells. 

The targeting of dyskerin is therefore expected to weaken both 
the production of ribosomes and the proper telomerase complex 
functioning impairing preferentially the growth of highly prolifer- 
ating cancer cells. Based on these facts we envisioned to specifically 
targeting dyskerin catalytic function by means of small molecule 
inhibitors in order to preferentially target cancer cells. To this end 
we have generated a structural model of the full-length human 
dyskerin based on known crystal structures of yeast [16] and 
screened, in silico, the Open National Cancer Institute database to 
select potential small molecule compounds interacting with the 
dyskerin catalytic site that was found particularly amenable for 
accommodating ligands [17,18]. 

The effects of the selected compounds were tested with a newly 
developed and optimized fast real-time PCR assay detecting the 
dyskerin pseudouridylation activity in vitro. Among all compounds, 
we identified four agents capable of significantly inhibit dyskerin 
activity. This result represents the proof of principle indicating that 
dyskerin catalytic activity can be modulated pharmacologically by 
means of small molecule inhibitors. 

Results and Discussion 

Computer-aided modeling and molecular design 

The three-dimensional structure of the human dyskerin 
structure has not been solved yet. In order to pursue a structure- 
based identification of dyskerin inhibitors, we created a homology 
model of human dyskerin by using a template crystal structure 
from saccharomyces cerevisiae (PDB ID: 3U28), which has a sequence 
identity of 73% with the human dyskerin sequence (Figure 1 and 
Figure SI in File SI) [16]. The next more similar structure was 
found to be the prokaryotic cbf5 of pyrococcus juriosus (PDBID: 
2EY4), which had a lower sequence alignment identity of 40% and 
therefore was not further considered [19]. The quality of the 
constructed model was assessed with the Procheck suite [20], and a 
molecular dynamics simulation that demonstrated the preserva- 
tion of the model protein folding with a maximum backbone 
RMSD fluctuation of 2.5A (Figure S2 in File SI). These results 
highlighted the good quality of the homology model obtained 
starting from the yeast structure (Figure 1). 

Encouraged by this, and with the aim to identify new small 
molecule inhibitors targeting the human dyskerin catalytic site, we 
used the three-dimensional coordinates of the obtained homology 
model to set up a structure-based virtual screening approach. By 
means of molecular docking techniques, we screened in silico the 
open NCI compounds dataset of the CoCoCo databases [18] and 
the first hundreds of hits were visually inspected against different 
criteria (see Materials and Methods). We selected 14 molecules 
(Figure S3 in File SI) that were acquired from the NCI for 
experimental testing. 

Effect of selected molecules on cellular 
pseudouridylation activity and viability 

To test the effect of the selected molecules on dyskerin 
pseudouridylation specific activity, we developed an in vitro real- 
time PCR-based assay (Figure 2A). A synthetic rRNA was used as 
a substrate and contained a small amount of the 28S rRNA 
sequence with uridines U4393/U4390, which are normally 
converted into a pseudouridine by dyskerin, and a 20 nucleotides 
long sequence at the 3 ' to specifically reverse- transcribe the rRNA 
substrate. The synthetic rRNA substrate was incubated with a 
nuclear cellular lysate to produce the pseudouridylated RNA 




Figure 1. Structure of the homology model of human dyskerin. 

The model is based on the template of saccharomyces cerevisiae 
crystallographic structure of the Cbf5-Nop10-Gar1 complex (PDB id: 
3U28). Important residues for pseudouridylation are highlighted to 
identify the catalytic region on the whole protein model. The present 
structure was obtained using MODELLER with the sequence alignment 
in Figure SI in File SI. The Figure was produced with Maestro. 
doi:1 0.1 371/journal.pone.01 01 971 .g001 

product. To detect substrate modification, the RNA product was 
chemically modified after lysate incubation using the [N- 
cyclohexyl-N , -P-(4-methyl-morpholinium)ethylcarbodiimide] 
(CMC) in order to generate N3-CMC-\|/ adducts. While CMC is 
known to bind U and G bases, these adducts are readily cleaved by 
weakly alkaline conditions. Conversely, cleavage of N3-CMC-\|/ 
requires strong alkaline conditions, and a combination of reactions 
with CMC and alkaline cleavage can be made specific to CMC 
addition to v|/. The obtained RNA product was subsequendy 
reverse-transcribed specifically. The CMC-\|/ complex efficiendy 
blocks reverse transcription creating a truncated PCR product, 
which is not amplified through the successive real-time PCR (as 
experimentally demonstrated in Figure S4 in File SI). As a 
consequence, the assay allows a semi-quantitative evaluation of 
dyskerin pseudouridylation activity by estimating the reverse- 
transcribed target products by real-time PCR. 

We first tested the efficiency of the in vitro assay (Figure 2B-D). 
In order to evaluate the capacity of our system to detect the 
presence of pseudouridines in selected positions, we tested our 
system using two different synthetic target RNA, one containing 
the uridines in positions 4393 and 4390 (U4393/4390) and one 
containing pseudouridines in the same positions (v|/4393/ 4390). To 
evaluate dyskerin pseudouridylation activity, both synthetic RNAs 
were incubated in a nuclear lysate at 30°C because this 
temperature represents the optimal condition for maximal 
enzymatic activity [21]. As an additional control we incubated 
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Synthetic rRNA substrate 28S4393/4390: 

(dGd)ATCCCGAACCTGAACCGGCCAAGCGuL)CAUAGCGACGCCCUAUAGUGAGUCGUAUUA(dTdC) 
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Figure 2. Pseudouridylation assay design and set up. (A) Schematic representation of the developed pseudouridylation real-time PCR-based 
assay. (8) Semi-quantitative evaluation of dyskerin pseudouridylation activity estimating the reverse-transcribed target products by real-time PCR. The 
same amount of two synthetic target RNAs, two different synthetic target RNA, one containing uridines in positions 4393 and 4390 and one 
containing pseudouridines in the same positions were subjected to nuclear lysate incubation, CMC modification, specific reverse transcription and 
real time PCR. Experiments were carried out at 4° and 30°C for 120 minutes. Histograms represent means and SDs from three independent 
experiments. *** = p<0.001). P value<0.05 is considered significant. (Q Semi-quantitative evaluation of dyskerin pseudouridylation activity in a 
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nuclear lysates from control (CTRL) and dyskerin depleted cells (shDKCI). Experiments were carried out at 30°C for 60 minutes. Histograms represent 
means and SDs from three independent experiments. *** = p<0.001 . P value<0.05 is considered significant. (D) Time course of the pseudouridylation 
reaction. The reactions were carried out at 30°C for 0, 15, 30, 60, 120, and 240 minutes and uridines U4393/U4390 were then verified by real-time PCR. 
Means and SDs from three independent experiments are represented. 
doi:1 0.1 371 /journal.pone.01 01 971 .g002 



the U4393/4390 synthetic rRNA substrate with the nuclear 
cellular lysate at 4°C, a condition in which the enzymatic reaction 
is down regulated. Our results showed that the rate of 
pseudouridylation of uridines U4393/U4390 at 30°C is significant 
increased (up to 12-fold) as respect to the same conditions at 4°C. 
On the other side, the comparison between the U4393/4390 and 
the i|/4393/4390 rRNAs allowed to demonstrate that a further 
increase in the level of pseudouridines can be measured by our 
system. (Figure 2B). Moreover, we evaluated the pseudouridylation 
activity by means of our assay in a nuclear lysate obtained from 
dyskerin depleted cells by DKC1 specific shRNA overexpression 
(shDKC 1) and in their relevant control [22] . Lysates from dyskerin 
depleted cells displayed a reduced pseudouridylation activity with 
respect to the control lysate (Figure 2C). These controls 
demonstrate that dyskerin pseudouridylation activity is measured 
by the in vitro assay developed. We then evaluated pseudouridyla- 
tion activity after 15, 30, 60, 120 and 240 minutes of the reaction 
incubation (Figure 2D). Results confirm that 120 minutes of 
incubation represent the best condition to appreciate the 
differences among the samples. To our knowledge, this is the first 
technology enabling the quantification of dyskerin catalytic activity 
without the use of radiolabelled tracers and the chromatographic 
separation of labeled nucleosides. 

We then tested the effect of selected molecules on dyskerin 
pseudouridylation activity using the assay described above 
(Figure 3A). In particular, we incubated each compound at 
100 (XM concentration in the reaction mixture for 120 minutes. 
The results showed that, among 14 tested compounds, compounds 
1, 5, 6 and 10 significantly reduced the specific pseudouridylation 
activity on 28S uridines U4393/U4390, although to different 
extents (Table 1 and Figure S3 in File SI). This result indicates, for 
the first time, the possibility to target the dyskerin catalytic site 
with small molecules. As a consequence, it is possible to postulate 
the use of small molecule inhibitors as a new approach to target 
dyskerin catalytic activity in tumors with an increased dyskerin 
expression. 

To test the biological effect of the four agents found, we 
evaluated their cytotoxicity using the Alamar Blue method. We 
performed the assay treating the MCF7 breast cancer cell line up 
to 96 hours. As shown in Figure 3B compounds 5, 6, and 10 had 
no cytotoxic effect in respect to control during the first 72 hours of 
treatment and the effect was detected only during the last 24 
hours. On the contrary, compound 1 had a significant cytotoxic 
activity starting from 24 h of treatment. The toxicity of compound 
1 was confirmed by the assessment of a dose/response curve 
(Figure S5 in File SI). 

By targeting the catalytic activity of dyskerin, we intended to 
affect preferentially uridine modification and not other non- 
catalytic functions of the protein such as TERC stabilization [9] . 
The obtained results indicated that, in the employed experimental 
conditions, compound 1 has no significant effect on TERC levels 
(Figure S6 in File SI). 

Moreover, we tested whether the effect of compound 1 on 
pseudouridylation activity is present not only on a nuclear lysate 
but on in vitro cultured cells. For this purpose we exposed MCF7 
cells to different concentrations of compound 1 and, adapting our 
RT-real time PGR approach, we measured on the endogenous 
modification of the same uridine pair (28S rRNA U4393/U4390) 



tested in the in cell free assay. Obtained results indicate that 
endogenous 28S rRNA U4393/U4390 modification is inhibited 
by compound 1 (Figure 3C). Importandy, the treatment with a 
compound displaying no effect on dyskerin pseudouridylation 
activity in the cell free assay (compound 9) did not induce any 
change on the modification of the same sites in cells (Figure S7 in 
File SI). 

Compound 1 is also known as pyrazofurin, a known inhibitor of 
the orotodine-5 '-monophosphate-decarboxylase (ODCase) in- 
volved in the de novo synthesis of pyrimidines nucleotides, DNA 
synthesis blocking and cellular replication [23,24]. As a conse- 
quence, it is possible to hypothesize that the cytotoxic activity 
observed was driven by the inhibition of both dyskerin and 
ODCase. The use of this compound has been reported in several 
clinical trials for ovarian carcinoma [25], sarcoma [26,27], 
colorectal carcinoma [28], acute myelogenous leukemia [29], 
breast cancer [30], lung cancer [31], melanoma [32] and others 
[33]. While these trials revealed a limited therapeutic potential of 
pyrazofurin, previous studies never considered dyskerin expres- 
sion. Thus, it is not possible to exclude the possibility that this 
compound may be highly active on a selected group of patients 
with dyskerin overexpression. In order to understand whether 
possible patterns could arise between overexpression profile of 
human dyskerin and bioactivity data, we envisioned highlighting 
possible correlations from reported data of compounds 1, 5, 6 and 
10 and the expression of DKC1 gene in human cancer cell lines 
(Table SI and S2 in File SI). Only for compound 1, pyrazofurin, 
we could retrieve a significant set of bioactivity data from the 
PubChem database. We found that, in human cancer cell lines 
overexpressing the DCK1 gene, there is a significant greater 
incidence of reported activity of pyrazofurin as respect to cell lines 
in which the DKC1 gene is underexpressed (p = 0.03) (Figure 3D). 
These results suggest that dyskerin expression may be an 
important factor in determining the activity of pyrazofurin. 

Compound structure and activity analysis 

The results presented in Figure 3 reveal for the first time the 
possibility for small molecules to interact with dyskerin catalytic 
activity. Thus, the identification of chemical scaffolds able to 
inhibit dyskerin function constitutes the first step towards the 
design of ad-hoc molecules with improved enzymatic and biological 
effects. In respect to this, and with the purpose to rationalize at a 
molecular level the bioactivity of compounds 1, 5, 6 and 10, we 
performed a thoughtful structural analysis of the molecules and a 
summary of the interactions in the catalytic site that are 
responsible to ligand - dyskerin binding. 

A more refined molecular docking technique was carried out to 
keep into account possible conformational changes in the active 
site residues, and to obtain a more accurate prediction of their 
binding modes. To this purpose an Induced Fit Docking (IFD) 
protocol was used with the human dyskerin model and the four 
active compounds. The IFD allows the movement of the residues 
in the binding site region leading to a more accurate prediction of 
binding modes. Compound 1 is structurally different from 
compounds 5, 6 and 10 (Table 1) and the IFD results demonstrate 
its ability to bind the catalytic site of dyskerin with different poses. 
However, due to the similarity of 1 to the uridine scaffold and the 
fact that, according to the X-ray structure of the box H/ ACA co- 
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Figure 3. Effect of selected molecules on cellular pseudouridylation activity and viability. (A) Bioactivity of the 14 selected compounds 
against pseudouridylation activity of human dyskerin in vitro. All the reactions were carried out at 30°C for 120 minutes with 100 uM of each 
compounds or DMSO. Histograms represent means and SEMs from three independent experiments. P values<0.05 are considered significant and are 
indicated with *(p<0.05), **(p<0.01) and ***(p<0.001). (fi) Cytotoxicity of compounds 1, 5, 6 and 10 on MCF7 cell line. 100 uM of compound 1, 5,6 
and 10 were added to cell and Alamar Blue assay performed each 24 hours for 4 days. Means and SDs from three independent experiments are 
represented. (Q Correlation between the expression profiles of DKC1 gene and reported bioactivity experiments of pyrazofurin (compound 1). Cells 
lines are divided into three groups according to dyskerin expression: overexpressed, underexpressed and not differentially expressed. Histograms 
represent the percentage of the cell lines in each group sensitive to the pyrazofurin cytotoxic activity. * = p<0.05. P value<0.05 is considered 
significant. The P is not indicated where correlations between the groups are not significant. (D) Effect of treatment with compound 1 on MCF7 
endogenous U4393/U4390 rRNA pseudouridylation. Cells were treated with 1 uM, 10 uM and 100 uM of compound 1 or DMSO for 24 hours. 
Histograms represent means and SEMs from three independent experiments. *** = p<0.001. P value<0.05 is considered significant. 
doi:1 0.1 371 /journal.pone.01 01 971 .g003 



crystallized with bound target-RNA, the pyrimidine moiety of 
uridine lies in the inner pseudouridylation catalytic site (PDB ID: 
3HAX), the correct binding pose of 1 is expected to be the one 
shown in Figure 4A [34]. Induced fit docking of compounds 5, 6 
and 10, unveil a similar superimposition of the uracil ring 
(Figure 4B-D) that performs similar interactions in the pseudour- 
idylation site. In particular, compound 1 presents two hydrogen 
bond interactions with the side chain of the catalytic D125, one 



with the pyrazolic hydroxyl group and the other with the amidic 
end. Conversely, compounds 5, 6 and 10 interact with the 
backbone nitrogen of D 125. The interactions with the side chain 
of the catalytic aspartate, D125, are particularly important to 
disrupt the formation of a salt bridge with the R227 (Figure 4), that 
previous structural analysis identified as the detector of the 
pseudouridylation conversion and subsequent substrate release 
[35,36]. For all active compounds, there are conserved hydrogen 
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Table 1. Identity of chemical compounds (chemical structures available in Figure S3 in File SI) and their relative pseudouridylation 
measurements. 



Relative pseudouridylation (Mean ± SEM) 



1 

5 
6 
10 



0.2304±0.1183 
0.8283±0.02631 
0.7101 ±0.0064 
0.7879±0.07261 



Compounds were tested at 100 uM concentration. 
doi:1 0.1 371 /journal.pone.01 01 971 .t001 



bond interactions with key backbone residues of T123, G223 and 
Y225 (Figure 4). Importantly, compound 1 can also interact with 
R248 and Y153 that, together with the different interaction 
pattern observed with D125, can explain the higher inhibitory 
activity in respect to 5, 6 and 10 (Figure 3A). Moreover, the IFD 
pose of compound 3 (Figure 4C) may also elucidate the slightly 
higher inhibitory activity in comparison to 5 and 10, since 
compound 3 is able to produce additional hydrogen bonds with 
Y153, G223 and Y225, by exploiting the ureic and hydroxyl 
groups direcdy linked to the uridine scaffold. 

It is worth to note that all these ligand-protein interactions 
involve polar residues with either charged or uncharged side 
chains. This fact highlights the polar nature of the binding pocket 



(Figure S8A in File SI) and also suggests the ideal region where a 
ligand feature should be located in order to interact optimally with 
the human dyskerin catalytic site. The ligand should mainly 
include hydrogen bonds donor and acceptor features (Figure S8B 
in File SI) [35]. The identification of the interaction patterns of the 
pseudouridylation catalytic site may be of practical use to define 
specific pharmacophore features useful for design of specific small- 
molecule modulators of human dyskerin. In fact, the comparative 
analysis of the available structures of related pseudouridine- 
synthase (e.g. PUS1 and PUS 10) showed that the catalytic site is 
extremely conserved at the structural level [37,38]. Therefore, it is 
possible that a certain degree of inhibitory capacity of the active 
compound 1 might be also present towards related enzymes 




Figure 4. Binding modes and interaction patterns of bioactive compounds on human dyskerin. The inhibitors are located in the 
pseudouridylation catalytic site of the human dyskerin homology model presented in Figure 1, The binding modes are the result of an Induced Fit 
docking protocol for each one of the four active compounds: (A) compound 1; (6) compound 5; (C) compound 6; (D) compound 10. Hydrogen bonds 
are highlighted in dotted green lines. Color scheme: white, hydrogen atoms; blue, nitrogen atoms; red, oxygen atoms; yellow, sulfur atoms; gray, 
dyskerin model carbon atoms; green, ligand carbon atoms. The Figure was produced with Maestro software. 
doi:1 0.1 371 /journal.pone.01 01 971 .g004 
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(Figure S9A in File SI). 28S rRNA U4393 and U4390 tested in 
our assays are known to be specifically modified by dyskerin 
through the binding of snoRNAs U68 and E3, respectively [39]. 
To further confirm that compound 1 effect on the in vitro assay is 
due to its inhibitory effect on dyskerin activity, we used a different 
synthetic substrate RNA depleted of the sequences for snoRNAs 
U68 and E3 recognition and for U4390, but conserving U4393. In 
this case, the results were not influenced by compound 1 (Figure 
S9B in File SI) indicating that the identified compound is certainly 
active towards dyskerin pseudouridylation activity. Nevertheless, 
further studies aiming at modifying the chemical structure of 
compound 1 are necessary to achieve target selectivity and, in this 
context, computational insights herein presented promise to 
constitute a useful tool to speed-up this development. 

Conclusions 

Our results indicate that it is possible to inhibit the dyskerin 
catalytic activity by means of low molecular weight compounds. In 
particular, we identified four molecules capable of significandy 
inhibit dyskerin pseudouridylation activity. The inhibition of the 
identified compounds is heterogeneous ranging from 18 to 77%. 
One of these compounds is pyrazofurin that was already tested as 
an anticancer drug. Because clinical trials with this agent did not 
take into account bioactivity data on human dyskerin, we cannot 
exclude that this compound may be readily for use in treating 
specific patients and/or cancer pathologies were dyskerin is 
overexpressed. 

Overall, the revelation that small molecule compounds can 
inhibit the catalytic function of dyskerin and the availability of a 
new fast assay to detect the dyskerin pseudouridylation activity 
in vitro herein described, will pave the way to a pharmacological 
targeting of ribosome biogenesis in cancer cells. In addition, the 
identification of the four chemical scaffolds described will serve as 
a starting point to optimize novel molecules with improved 
pharmacological properties and biological profiles against human 
dyskerin. 

Methods 

Homology modeling 

A homology model was created starting from the human 
dyskerin protein sequence (Uniprot code: 060832), by excluding 
nuclear and nucleolar localization sequences, which are not 
related to the target pseudouridylation catalytic site [40]. The 
resulting primary sequence (Figure SI in File SI) was submitted to 
BLAST by searching on the Protein Data Bank (PDB) for the 
possible existence of three-dimensional homologous proteins of the 
human dyskerin [41]. The best homology score was found for the 
Cbf5-Nopl0-Garl complex from saccharomyces cerevisiae (PDB ID: 
3U28), with a sequence alignment identity of the 73% [16]. The 
pair-wise sequence alignment was performed with ClustalW and 
manually inspected (Figure SI in File SI) [42]. 

We created the homology model of human dyskerin using the 
yeast template with MODELLER software (Figure 1 and Figure 
SI in File SI) [43]. The model was refined and submitted to a 
quality check with the Procheck suite. Subsequendy, a molecular 
dynamics simulation was performed in order to assess the quality 
and the folding of the constructed human dyskerin model. To this 
purpose the Desmond software (v3. 1.5 1.1) was used in construct- 
ing a SPC water simulation box with standard parameters and to 
run the molecular simulation for 10 ns. 



Virtual screening 

With the aim to discover new small molecule inhibitors 
targeting the human dyskerin catalytic site, the docking software 
Glide (version 5.7), was used to perform a high-throughput virtual 
screening of the open NCI (http://dtp.cancer.gov) compound 
dataset of the CoCoCo databases [18]. A total of around 260.000 
molecules were screened by centering a 20 A docking grid in the 
D125 residue, which has been described as the amino acid 
responsible of the pseudouridylation function of dyskerin 
[19,35,44]. Docking results were ranked based on the Glide score 
and the first 500 hit compounds were visually inspected taking into 
account several structural and physic-chemical rules, such as the 
qualitative evaluation of ligand-protein interactions within the 
active site, probability of suggested protonation and tautomeric 
states, stereochemistry complexity, compound availability and 
chemical diversity. Based on these criteria, 14 molecules (Table 1 
and Figure S3 in File SI), all listed on the top of the 200 scoring 
hit, were acquired from the NCI for experimental testing. 

Induced fit docking 

The Induced fit docking protocol present in Schrodinger Suite 
2012 is a procedure that takes into account the protein flexibility in 
a protein-ligand docking [45,46] . To achieve more refined docking 
results this protocol includes four major steps: an initial softened- 
potential docking, a sampling of the protein in the binding region, 
and finally a redocking of the ligand into low energy induced-fit 
structures resulted from the residues sampling followed by a 
scoring that reports the docking energy (GlideScore), and the 
receptor strain and solvation terms (Prime energy). The Induced fit 
docking protocol was applied to the human dyskerin model 
previously described and each one of the active compounds: 1, 5, 
6 and 10. 

Pseudouridylation assay 

Synthetic rRNA substrate 28S U4393/U4390 design. The 

synthetic rRNA substrate 28S U4393/U4390 was designed to 
contain a 22 nucleotides long sequence of 28S rRNA with uridines 
U4393/U4390 and a 20 nucleotides long sequence at the 3' to 
specifically reverse-transcribe the rRNA substrate. The complete 
sequence is here reported: (dGd)ATCCCGAACCT- 
GAAC C GGC CAAGC GUUC AUAGC GAC GC C CUAUAGU- 
GAGUCGUAUUA(dTdC) 

A similar synthetic sequence containing pseudouridines in 
positions corresponding to uridines 4393/4390 in 28S RNA was 
used as a control. The RNA oligonucleotides were purchased from 
IDT DNA. 

Nuclear lysate preparation. Nuclear lysate was prepared 
from four subconfluent T175 flasks of MCF7 cells. The cells were 
scraped off, washed twice with 10 ml of cold phosphate-buffered 
saline (PBS), and swollen for 1 5 min in 2 times the pellet volume of 
a hypoosmotic homogenization buffer (10 mM HEPES pH 7.5; 
1.5 mM MgC12; 10 mM KC1 and 1 mM DTT). The cells were 
disrupted with 12 to 15 strokes of a tight-fitting pestle in a Dounce 
homogenizer. After a 10-min centrifugation at 1,000 g, the pellets 
(crude nuclei) were suspended in 0.5 ml of reconstitution buffer 
(10 mM Tris-HCl pH 7.5; 3,3 mM MgC12; 0,25 M sucrose and 
1 mM DTT). The viscous solution was sonicated 6-8 times for 
15 sec followed by 15-sec cooling periods in between until a great 
reduction of the viscosity occurred. The nuclear lysate was frozen 
in aliquots and stored at — 80C. 

Pseudouridylation reaction. Pseudouridylation reaction 
was carried out as previously described (20) with some modifica- 
tions. Particularly, the reaction was carried out in a final volume of 
200 (Xl in order to obtain a sufficient amount of RNA to be used in 
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the consequent passages. A 200 JLll reaction contained 80 JLll of 
MCF7 nuclear extract, 40 Ul of 5X Reconstitution buffer 
(100 mM HEPES-KOH pH 7.5; 600 mil KC1; 10 raM MgCl 2 ), 
8 pg of tRNA, 1 60 U of RNasin (Promega) and approximately 80 
fmol of rRNA substrate. Each compound was diluted in DMSO 
(Sigma-Aldrich) and 100 (iM of each of them or DMSO (in the 
control samples) were added to the mixture. The mixture was 
incubated at 30° and 4°C for 120 minutes unless differendy 
indicated. The reaction was terminated on ice and RNA extracted 
with Tri- Reagent (Ambion) following the manufacturer's proto- 
col. 

CMC-modification of rRNA. CMC modification was car- 
ried out as previously described with some modifications [47] . To 
2 Ug of RNA were added 80 ul of BEU Buffer (7 M Urea, 4 mM 
EDTA, 50 mM Bicine-pH 8.5, final buffer pH = 8.9-9.0) and 
20 Ul of 1 M CMCT (Sigma-Aldrich) in BEU buffer freshly 
prepared. The samples were incubated for 20 minutes at 37°C 
and then precipitated by adding 2 (Lll of Pellet Paint Co-Precipitant 
(Merck), 50 ul of 3 M sodium acetate pH 5.5 and 600 ul of 
ethanol. The pellets were washed twice with 200 ul of 70% 
ethanol and then resuspended in 50 |il sodium carbonate buffer, 
pH 10.4 (50 mM sodium carbonate and 2 mM EDTA). The 
samples were then incubated at 37°C for 4 h. At the end of the 
reaction RNA was precipitated by adding 2 Ul of Pellet Paint Co- 
Precipitant, 6 ul of 3 M sodium acetate, and 110 ul of ethanol. 
The pellet was then washed with 70% ethanol, resuspended in 
20 Ul of water and finally quantified by a nanodrop. 

Reverse transcription reaction. The reverse transcription 
reaction was carried out using a specific reverse primer targeting 
the 20 nucleotides long sequence in the 3' region of the synthetic 
rRNA substrate. The sequence of the specific reverse target is: 5'- 
TAATACGACTCACTATAGGG-3 ' . 

The reaction was made in two steps in a final 20 (J.1 mixture. 
During the first step 300 ng of RNA, 1 uM specific reverse primer 
(IDT DNA) and 40 U of RNasin Plus RNase Inhibitor (Promega) 
were mixed and incubate at 50° for 5 minutes in a thermocycler to 
stimulate the annealing of the specific reverse primer. During the 
second step, M-MLV Reverse Transcriptase Reaction Buffer 
(Promega), 200 U of M-MLV retrotranscriptase (Promega) and 
10 UM dNTPs (Promega) were added to mixture and inserted in 
the thermocycler. The conditions of the reverse transcription were: 
5 minutes at 4°C, 60 minutes at 37°C, 15 minutes at 70°C and 
4°C forever. 

Real time PCR reaction. Real-time PCR analysis was 
performed in a Gene Amp 7000 Sequence Detection System 
(Applied Biosystems) using the Syber-green (Applied Biosystems) 
approach. The set of primers (IDT DNA) used were: primerF 5'- 
GAACCTGAACCGGCCAAG-3' and primerR 5'-CTA- 
TAGGGCGTCGCTATGA-3 ' . This set of primer permits the 
amplification of the specific retro transcribed RNA. For each 
sample three replicates were analyzed. The amount of 28S 
U4393/U4390 uridines pseudouridylated was calculated by the 
fold change 2 A DCt were DCt is the difference between the Ct 
medium of the samples treated with each molecule and the Ct 
medium of the positive control (the sample without the addition of 
inhibitors at 30°C). Real Time PCR for human TERC mRNA 
expression analysis was performed as previously described [48,49] , 
using human U6 snRNA, as endogenous control. Primers for U6 
snRNA expression analysis were the following: primerF 5'-GCT 
GGC TTC GGC AGC ACA TAT AC-3', primerR 5'- 
TATCGAACGCTTCACGAATTTGC-3'. The final results were 
determined by the 2 A — AACt method. 



Cell cultures and viability assay 

Human breast cancer derived cell lines MCF7 (obtained from 
ATCC) were cultured in a monolayer at 37°C in a humidified 
atmosphere containing 5% C0 2 . MCF7 were grown in RPMI 
1640 (Sigma-Aldrich) supplemented with 10% of fetal bovine 
serum (FBS, Sigma-Aldrich)), 2 mM L-glutamine (Sigma-Al- 
drich)), 100 U/ml penicillin, and 100 mg/ml streptomycin (Sig- 
ma-Aldrich). Cells were treated with 100 uM of each compound 
or DMSO for 24, 48, 72 and 96 hour. For the dose/response 
curve cells were treated with 5, 10, 50 and 100 uM of compound 1 
for 24, 48, 72 and 96 hours. We used the Alamar Blue assay 
(Invitrogen) to test the vitality of the cells after the treatment 
following the manufacturer's instruction. Stable MCF7 (shDKCl), 
and their relevant empty vector controls were generated as 
described previously [22]. 

Pseudouridylation assay on endogenous 28S rRNA 

MCF7 cells were treated with 1 uM, 10 uM and 100 uM of 
compound 1, 9 or DMSO for 24 hours. Two Ug of RNA was 
modified with CMC as described above. The reverse transcription 
reaction was carried out using an oligonucleotide directed to a 
sequence adjacent to the same 28S RNA sequence included in the 
synthetic rRNA substrate used in the in vitro assay: 5'- 
ATTATGCTGAGTGATATCCCATCGAAGGAT- 
CAAAAAGCGA-3'. The oligonucletide contained a linker 
sequence in its 5 ' allowing the specific amplification of the reverse 
transcribed templates. The reaction was made in two steps in a 
final 20 ul mixture. During the first step 100 ng of RNA, 0.2 uM 
specific reverse primer and 40 U of RNasin Plus RNase Inhibitor 
(Promega) were mixed and incubate at 50°C for 10 minutes to 
allow the annealing of the specific reverse primer. During the 
second step, M-MLV Reverse Transcriptase Reaction Buffer 
(Promega), 200 U of M-MLV retrotranscriptase (Promega) and 
2 uM dNTPs (Promega) were added to mixture. Reverse 
transcription and Real-time PCR analysis was performed as 
described above, using the following PCR primers: primer F 5'- 
GTGTC AG AAAAGTTAC C AC A- 3 ' ; primer R: 5'-AT- 
TATGCTGAGTGATATCCC-3'. 

Statistical analysis 

Mann-Whitney U or Fischer's exact tests were used, when 
appropriate, for the comparisons among groups. Values for p less 
than 0.05 were regarded as statistically significant. 

Expression profile analysis 

Activity and inactivity data of compound 1 were retrieved from 
the PubChem database (accessed Feb, 21 2013 http://pubchem. 
ncbi.nlm.nih.gov; CID 285701). The number of published studies 
where DKC1 gene is differentially expressed compared to the 
gene's overall mean expression level in the study was retrieved 
from the Gene Expression Atlas database (accessed Feb, 21 2013 
http://www.ebi.ac.uk; DKC1 homo sapiens). 

Supporting Information 

File SI Supporting Information. Methods SI, Pseudour- 
idylation assay on synthetic rRNA lacking E3/U68 binding 
sequences. Figure SI, Sequence alignment of human dyskerin 
and yeast cbf5 primary structures. Figure S2, Backbone RMSD 
fluctuation of the human dyskerin model presented and used for 
structure based virtual screening during 10 ns of molecular 
dynamics simulation. The higher fluctuation between 7-8.5 ns is 
due to the movement of the loop covering the pseudouridylation 
site. The protein remains folded for the complete simulation time 
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demonstrating the good quality of the assembled human dyskerin 
model. Figure S3, Chemical structures of tested compounds for 
inhibition of pseudouridylation catalytic activity of human 
dyskerin. Figure S4, Real time RT-PCR assay used in the 
pseudouridylation assay selectively amplifies full length products 
(5' TAA TAC GAC TCT CTA TAG GGC GTC GCT ATG 
AAC GCT TGG CCG GTT CAG GTT CGG GAT 3'), while 
predicted truncated products (5' TAA TAC GAC TCT CTA 
TAG GGC GTC GCT 3') originating from pseudouridylated 
templates are not amplified. Figure S5, Dose-Response curve of 
compound 1 on MCF7 cells. The experiment was performed 
treating the cells with DMSO (CTRL) or 5, 10, 50 and 100 uM of 
compound 1 for up to 96 hours. Alamar Blue assay was performed 
each 24 hours. Means and SDs from three independent 
experiments are represented. Figure S6, Effect of compound 1 
on telomerase RNA component (TERC) expression. The 
experiment was performed treating the nuclear lysate with DMSO 
(CTRL) or 100 (XM of compound 1 for 120 minutes in the 
reaction mixture. Histograms represent means and SDs from three 
independent experiments. The final results were determined by 
the 2 A — AACt method. Differences between the groups are not 
significant. Figure S7, Effect of treatment with compound 9 on 
MCF7 endogenous U4393/U4390 rRNA pseudouridylation. 
Cells were treated with 1 [iM, 10 |0M and 100 |0,M of compound 
9 or DMSO for 24 hours. The pseudouridylation reaction was 
carried out at 30°C. Histograms represent means and SEMs from 
three independent experiments. Differences between the groups 
are not significant. Figure S8, (A) Contour maps generated with 
the software SiteMap (http://www.schrodinger.com), thereby 
presenting hydrophobic (yellow regions), donor (blue regions) 
and acceptor (red regions) potentials. Contour maps represent the 
ideal region of the space where a corresponding ligand feature 
should be located in order to interact optimally with the human 
dyskerin catalytic site. (B) Molecular surface of the human dyskerin 
coloured by electrostatic potential. Blue regions represent 
positively charged residues while red regions negatively charged 
residues. Figure S9, (A) Comparative analysis of the available 

References 

1 . Montanaro L, Trerc D, Dcrcnzini M (2008) Nucleolus, ribosomcs, and cancer. 
Am J Pathol 173: 301-310. Available: http://www.pubmcdccntral. nih.gov/ 
articlercndcr.fcgi?artid — 2475768&tool — pmccntrcz&rcndcrtypc — abstract. Ac- 
cessed 10 February 2013. 

2. Montanaro L, Trerc D, Dcrcnzini M (2012) Changes in ribosome biogenesis 
may induce cancer by down-rcgulating the cell tumor suppressor potential. 
Biochim Biophys Acta 1825: 101-110. Available: http://www.ncbi.nlm.nih. 
gov/pubmed/22079382. Accessed 10 February 2013. 

3. Wright WE, Piatyszck MA, Raincy WE, Byrd W, Shay JW (1996) Telomerase 
activity in human gcrmlinc and embryonic tissues and cells. Dev Genet 18: 173— 
179. Available: http://www.ncbi. nlm.nih.gov/pubmed/8934879. Accessed 10 
February 2013. 

4. Collins K, Mitchell JR (2002) T elomerase in the human organism. Oncogene 
21: 564-579. Available: http://www.ncbi.nlm.nih.gov/pubmed/11850781. 
Accessed 10 February 2013. 

5. Drygin D, Rice WG, Grummt I (2010) The RNA polymerase I transcription 
machinery: an emerging target for the treatment of cancer. Annu Rev 
Pharmacol Toxicol 50: 131-156. Available: http://www.ncbi.nlm.nih.gov/ 
pubmcd/20055700. Accessed 10 February 2013. 

6. Drygin D, Siddiqui-Jain A, O'Brien S, Schwaebc M, Lin A, ct al. (2009) 
Anticancer activity of CX-3543: a direct inhibitor of rRNA biogenesis. Cancer 
Res 69: 7653—7661. Available: http://www.ncbi.nlm.nih.gov/pubmed/ 
19738048. Accessed 10 February 2013. 

7. Harley CB (2008) Telomerase and cancer therapeutics. Nat Rev Cancer 8: 167- 
179. Available: http://www.ncbi.nlm.nih.gov/pubmcd/18256617. Accessed 10 
February 2013. 

8. Shay JW, Wright WE (2002) Telomerase: a target for cancer therapeutics. 
Cancer Cell 2: 257-265. Available: http://www.ncbi. nlm.nih.gov/pubmed/ 
12398889. Accessed 30 January 2013. 

9. Montanaro L (2010) Dyskerin and cancer: more than telomerase. The defect in 
mRNA translation helps in explaining how a proliferative defect leads to cancer. 



structures of related pseudouridine-synthase PUS1 and PUS 10. 
On the left column, compound 1 complexes were minimized into 
the active site of PUS 1 and PUS 10 starting from the binding pose 
obtained by overlaying these structures with the human Dyskerin - 
compound 1 complex (Yellow lines indicate hydrogen bonds); On 
the right column, a ligand interaction diagram for each complex is 
presented (Grey highlight - solvent exposed; blue arrow - side 
chain hydrogen bond; blue dotted arrow - backbone hydrogen 
bond; red ball - negative charged residue; blue ball - positive 
charged residue; green ball - hydrophobic residue; light blue ball - 
polar residue). (B) Specific inhibitory effect of compound 1 on 
dyskerin activity. A different synthetic substrate RNA depleted of 
the sequences for snoRNAs U68 and E3 recognition and for 
U4390, but conserving U4393 was used in the in vitro assay 
developed. The experiment was performed treating the nuclear 
lysate with DMSO (CTRL) or 100 uM of compound 1 for 
1 20 minutes in the reaction mixture. Histograms represent means 
and SEMs from three independent experiments. Correlations 
between the groups are not significant. Table SI, List of NCI 
human tumor cell lines with confirmed bioactivity of compound 1 
and their relative expression profiles of human dyskerin classified 
as overexpressed, underexpressed or not differentially expressed. 
Data were collected as described in the materials and methods 
paragraph. Table S2, List of NCI human tumor cell lines with 
confirmed inactivity of compound 1 and their relative expression 
profiles of human dyskerin classified as overexpressed, under- 
expressed or not differentially expressed. Data were collected as 
described in the materials and methods paragraph. 
(PDF) 

Author Contributions 

Conceived and designed the experiments: ADR LM. Performed the 
experiments: LR AJMB CO. Analyzed the data: LR AJMB CO ADR LM. 
Contributed reagents/materials/analysis tools: ADR LM. Wrote the 
paper: LR AJMB CO ADR LM. 



J Pathol 222: 345—349. Available: http://www.ncbi. nlm.nih.gov/pubmcd/ 
20925138. Accessed 30 March 2012. 

10. Montanaro L, Galicnni M, Bertoni S, Rocchi L, Sansone P, et al. (2010) Novel 
dyskcrin-mcdiatcd mechanism of p53 inactivation through defective mRNA 
translation. Cancer Res 70: 4767-4777. Available: http://www.ncbi.nlm.nih. 
gov/pubmed/20501855. Accessed 2 February 2013. 

11. Sieron P, Hadcr C, Hatina J, Engers R, Wlazlinski a, et al. (2009) DKC1 
overexpression associated with prostate cancer progression. Br J Cancer 101: 
1410-1416. Available: http:/ /www. pubmedcentral.nih.gov/ articlerendcr. 
fcgi?artid - 2768451&tool — pmccntrcz&rcndcrtypc = abstract. Accessed 10 
February 2013. 

12. Bellodi C, Krasnykh O, Haynes N, Theodoropoulou M, Peng G, et al. (2010) 
Loss of function of the tumor suppressor DKC 1 perturbs p27 translation control 
and contributes to pituitary tumor igene sis. Cancer Res 70: 6026-6035. 
Available : http: / / www. pubmedcentral.nih.gov / articlerender. 
fcgiPartid = 2913864&tool — pm centre z&rendcrtype = abstract. Accessed 23 Jan- 
uary 2014. 

13. Alawi F, Lin P, Ziober B, Patel R (201 1) Correlation of dyskerin expression with 
active proliferation independent of telomerase. Head Neck 33: 1041-1051. 
Available : http : / / www.pubmedcentral.nih.gov / articlerender. 
fcgi?artid = 31 16033&tool — pmccntrcz&rcndcrtypc = abstract. Accessed 10 
February 2013. 

14. Turano M, Angrisani A, De Rosa M, Izzo P, Furia M (2008) Real-time PCR 
quantification of human DKC1 expression in colorectal cancer. Acta Oncol 47: 
1598-1599. Available: http://www.ncbi.nlm.nih.gov/pubmed/18607840. Ac- 
cessed 10 February 2013. 

15. Liu B, Zhang J, Huang C, Liu H (2012) Dyskerin overexpression in human 
hepatocellular carcinoma is associated with advanced clinical stage and poor 
patient prognosis. PLoS One 7: c43147. Available: http://www.pubmedccntral. 
nih.gov / articlerender. 

fcgi?artid = 3418259&tool — pmccntrcz&rcndcrtypc = abstract. Accessed 10 
February 2013. 



PLOS ONE | www.plosone.org 



9 



July 2014 | Volume 9 | Issue 7 | e101971 



Targeting Dyskerin and Ribosome Biogenesis 



16. Li S, DuanJ, Li D, Yang B, Dong M, ct al. (201 1) Rcconstitution and structural 
analysis of the yeast box H/ACA RNA-guidcd pscudouridinc synthase. Genes 
Dev 25: 2409-2421. Available: http://www.pubmedcentral.nih.gov/ 
articlcrcndcr.fcgiPartid = 3222906&tool = pmcentrez&rendertype = abstract. Ac- 
cessed 10 February 2013. 

17. Moura Barbosa AJ, Del Rio A (2012) Freely accessible databases of commercial 
compounds for high- throughput virtual screenings. Curr Top Med Chem 12: 
866-877. Available: http://www.ncbi.nlm.nih.gov/pubmed/22352914. Ac- 
cessed 20 March 2013. 

18. Del Rio A, Barbosa AJM, Gaporuscio F, Mangiatordi GF (2010) GoCoCo: a free 
suite of multiconfbrmational chemical databases for high-throughput virtual 
screening purposes. Mol Biosyst 6:21 22-2 1 28. Available: http: / / www.ncbi.nlm. 
nih.gov/pubmcd/20694263. Accessed 14 March 2012. 

19. Rashid R, Liang B, Baker DL, Youssef O, He Y, et al. (2006) Crystal structure of 
a Cbf5-Nopl0-Garl complex and implications in RNA-guidcd pscudouridyla- 
tion and dyskeratosis congenita. Mol Cell 21: 249-260. Available: http://www. 
ncbi.nlm.nih.gov/pubmed/16427014. Accessed 30 March 2012. 

20. Laskowski RA, MacArthur MW, Moss DS, Thornton JM (1993) PROGHECK: 
a program to check the stereochemical quality of protein structures. J Appl 
Crystallogr 26: 283-291. Available: citculikc-articlc-id: 4778677. 

21. Wang C, Query CC, Meier UT (2002) Immunopurificd Small Nucleolar 
Ribonuelcoprotcin Particles Pscudouridylatc rRNA Independently of Their 
Association with Phosphorylated Noppl40. Mol Cell Biol 22: 8457-8466. 
Available: http://mcb.asm.org/egi/doi/10.1128/MCB.22.24.8457-8466.2002. 
Accessed 10 February 2013. 

22. Rocchi L, Pacilli A, Sethi R, Penzo M, Schneider RJ, ct al. (2013) Dyskerin 
depletion increases VEGF mRNA internal ribosome entry site-mediated 
translation. Nucleic Acids Res: 1—11. Available: http://www.ncbi.nlm.nih. 
gov/pubmcd/23821664. Accessed 9 August 2013. 

23. Cadman EC, Dix DE, Handschumacher RE (1978) Clinical, biological, and 
biochemical effect of pyrazofiirin. Cancer Res 38: 682-688. Available: http:// 
www.ncbi.nlm.nih.gov/pubmcd/272228. Accessed 10 February 2013. 

24. Christopherson RI, Lyons SD, Wilson PK (2002) Inhibitors of de novo 
nucleotide biosynthesis as drugs. Acc Chem Res 35: 96 1-97 1 . Available: 
http://www.ncbi.nlm.nih.gov/pubmcd/12437321. Accessed 10 February 2013. 

25. EdmonsonJH, Decker DC, Malkasian CD, Webb MJ (n.d.) Concomitant phase 
II studies of pyrazofurin and razoxane in alkylating agent-resistant cases of 
epithelial ovarian carcinoma. Cancer Treat Rep 65: 1127—1129. Available: 
http://www.ncbi.nlm.nih.gov/pubmed/7296557. Accessed 10 February 2013. 

26. Cormier WJ, Hahn RG, Edmonson JH, Eagan RT (n.d.) Phase II study in 
advanced sarcoma: randomized trial of pyrazofurin versus combination 
cyclophosphamide, doxorubicin, and eis-dichlorodiammineplatinum(II) (CAP). 
Cancer Treat Rep 64: 655-658. Available: http:/ / www.ncbi.nlm. nih.gov/ 
pubmcd/7191779. Accessed 10 February 2013. 

27. Gralla RJ, Sordillo PP, Magill GB (1978) Phase II evaluation of pyrazofurin in 
patients with metastatic sarcoma. Cancer Treat Rep 62: 1573-1574. Available: 
http://www.ncbi. nlm.nih.gov/pubmed/361228. Accessed 10 February 2013. 

28. Carroll DS, Kemeny NE, Gralla RJ (1979) Phase II evaluation of pyrazofurin in 
patients with advanced colorectal carcinoma. Cancer Treat Rep 63: 139-140. 
Available: http://www.ncbi.nlm.nih.gov/pubmed/369686. Accessed 10 Febru- 
ary 2013. 

29. Voglcr WR, Trulock PD (1978) Phase I study of pyrazofurin in refractory acute 
myelogenous leukemia. Cancer Treat Rep 62: 1569-1571. Available: http:// 
www.ncbi.nlm.nih.gov/pubmed/152146. Accessed 10 February 2013. 

30. Nichols WC, Kvols LK, Ingle JN, EdmonsonJH, Ahmann DL, ct al. (1978) 
Phase II study of triazinate and pyrazofurin in patients with advanced breast 
cancer previously exposed to cytotoxic chemotherapy. Cancer Treat Rep 62: 
837-839. Available: http://www.ncbi.nlm.nih.gov/pubmed/657166. Accessed 
10 February 2013. 

31. Gralla RJ, Currie VE, Wittes RE, Golbey RB, Young CW (1978) Phase II 
evaluation of pyrazofurin in patients with carcinoma of the lung. Cancer Treat 
Rep 62: 451-452. Available: http://www.ncbi.nlm.nih.gov/pubmed/348315. 
Accessed 10 February 2013. 

32. Budman D, Curric V, Wittes R (1977) Phase II trial of pyrazofurin in malignant 
melanoma. Cancer Treat Rep 61: 1733-1734. Available: http://www.ncbi.nlm. 
mh.gov/pubmed/340040. Accessed 10 February 2013. 

33. Cummings FJ, Stollcr RG, Calabresi P (1979) Clinical trial of weekly 
pyrazofurin. Cancer Treat Rep 63: 1363—1365. Available: http://www.ncbi. 
nlm.nih.gov/pubmed/157807. Accessed 10 February 2013. 



34. DuanJ, Li L, Lu J, Wang W, Ye K (2009) Structural mechanism of substrate 
RNA recruitment in H/ACA RNA-guidcd pscudouridinc synthase. Mol Cell 34: 
427-439. Available: http://www.ncbi.nlm.nih.gov/pubmed/19481523. Ac- 
cessed 11 February 2013. 

35. Hamma T, Ferre-D'Amarc AR (2010) The box H/AGA ribonucleoprotein 
complex: interplay of RNA and protein structures in post-transcriptional RNA 
modification. J Biol Chem 285: 805-809. Available: http: //www. 
pubmedccntral.nih.gov/articlerendcr. 
fcgi?artid — 2801 282&tool — pmcentrez&rendertype — abstract. Accessed 12 
February 2013. 

36. Hoang C, Hamilton CS, Mueller EG, Ferre-D'Amarc AR (2005) Precursor 
complex structure of pscudouridinc synthase TruB suggests coupling of active 
site perturbations to an RNA-sequestering peripheral protein domain. Protein 
Sci 14: 220 1-2206. Available: http:/ /www. pubmedeentral.nih.gov/ 
article render. fcgi?artid = 2279332&tool = pmcentrez&rendertype - abstract. Ac- 
cessed 20 February 2013. 

37. Czudnochowski N, Wang AL, Finer-Moore J, Stroud RM (2013) In human 
pseudouridinc synthase 1 (hPusl), a C -terminal helical insert blocks tRNA from 
binding in the same orientation as in the Pusl bacterial homologue TruA, 
consistent with their different target selectivities. J Mol Biol 425: 3875-3887. 
Available: http://www.ncbi.nlm.nih.gov/pubmed/23707380. Accessed 28Jan- 
uary 2014. 

38. Hamma T, Ferre-D'Amarc AR (2006) Pscudouridinc synthases. Chem Biol 13: 
1125-1135. Available: http://www.ncbi.nlm.nih.gov/pubmed/17113994. Ac- 
cessed 28 January 2014. 

39. Ganot P, Bortolin ML, Kiss T (1997) Site-specific pscudouridinc formation in 
preribosomal RNA is guided by small nucleolar RNAs. Cell 89: 799-809. 
Available: http://www.nebi.nlm.nih.gov/pubmed/9182768. Accessed 23Janu- 
ary 2014. 

40. Consortium TU (2012) Reorganizing the protein space at the Universal Protein 
Resource (UniProt). Nucleic Acids Res 40: D71-5. Available: http://www. 
pubmedccntral.nih.gov/articlerendcr. 
fcgi?artid- 32451 20&tool = pmcentrez&rendertype = abstract. Accessed 9 Feb- 
ruary 2013. 

41. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, et al. (2000) The 
Protein Data Bank. Nucleic Acids Res 28: 235-242. Available: http://www. 
pubmedccntral.nih.gov/articlerendcr. 
fcgi?artid = 102472&tool = pmcentrez&rendertype = abstract. 

42. Goujon M, McWilliam H, Li W, Valentin F, Squizzato S, et al. (2010) A new 
bioinformatics analysis tools framework at EMBL-EBI. Nucleic Acids Res 38: 
W695-9. Available: http:/ /www. pubmedccntral.nih.gov/ articlerender. 
fcgiPartid = 2896090&tool — pmcentrez&rendertype = abstract. Accessed 20 
February 2013. 

43. Sali A, Blundell TL (1993) Comparative protein modelling by satisfaction of 
spatial restraints. J Mol Biol 234: 779—815. Available: http:/ / www.ncbi.nlm.nih. 
gov/pubmed/8254673. Accessed 11 February 2013. 

44. Zebarjadian Y, King T, Fournier MJ, Clarke L, Carbon J (1999) Point 
mutations in yeast CBF5 can abolish in vivo pscudouridylation of rRNA. Mol 
Cell Biol 19: 7461-7472. Available: http://www.pubmedcentral.nih.gov/ 
articlercndcr.fegi?artid — 84741 &tool — pmcentrez&rendertype — abstract. Ac- 
cessed 20 February 2013. 

45. Sherman W, Day T, Jacobson MP, Fricsncr RA, Farid R (2006) Novel 
procedure for modeling ligand/rcecptor induced fit effects. J Med Chem 49: 
534-553. Available: http://www.ncbi.nlm.nih.gov/pubmed/16420040. Ac- 
cessed 11 February 2013. 

46. Sherman W, Beard HS, Farid R (2006) Use of an induced fit receptor structure 
in virtual screening. Chem Biol Drug Des 67: 83-84. Available: http://www. 
ncbi.nlm.nih.gov/pubmcd/16492153. Accessed 11 February 2013. 

47. Ofengand J, Del Campo M, Kaya Y (2001) Mapping pseudouridines in RNA 
molecules. Methods 25: 365—373. Available: http://www.ncbi.nlm.nih.gov/ 
pubmed/1 1860291. Accessed 20 February 2013. 

48. Montanaro L, Brigotti M, ClohcssyJ, Barbieri S, Ceecarelli C, ct al. (2006) 
Dyskerin expression influences the level of ribosomal RNA pseudo-uridylation 
and telomcrase RNA component in human breast cancer. J Pathol 210: 10-18. 
Available: http://www.ncbi.nlm.nih.gov/pubmed/16841302. Accessed 23Jan- 
uary2014. 

49. Montanaro L, Calienni M, Ceecarelli C, Santini D, Taffurclli M, ct al. (2008) 
Relationship between dyskerin expression and telomcrase activity in human 
breast cancer. Cell Oncol 30: 483—490. Available: http://www.ncbi.nlm.nih. 
gov/pubmed/ 18936525. Accessed 23 January 2014. 



PLOS ONE | www.plosone.org 



10 



July 2014 | Volume 9 | Issue 7 | e101971 



